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Abstract 



The Ab differential production cross section and the cross-section ratio cr(Ab)/c(Ab) 
are measured as functions of transverse momentum pj b and rapidity |y Ab | in pp colli- 
sions at i/s = 7 TeV using data collected by the CMS experiment at the LHC The mea- 
surements are based on Ab decays reconstructed in the exclusive final state ]/ipA, with 
the subsequent decays ]/tp — > and A — > Tpn, using a data sample corresponding 

to an integrated luminosity of 1.9fb _1 . The product c(Ab) x £>(Ab — > J/tpA) versus 
Pj h falls faster than that of b mesons. The measured value of c(Ab) x £>(Ab — > ]/ipA) 
for pj b > lOGeV and \y A *\ < 2.0 is 1.16 ± 0.06 ± 0.12 nb, and the integrated 
cr(Ab)/c(Ab) ratio is 1.02 ± 0.07 ± 0.09, where the uncertainties are statistical and 
systematic, respectively. 
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1 Introduction 

Cross sections for b-quark production in high-energy hadronic collisions have been measured 
at pp colliders at center-of-mass energies from 630 GeV HI to 1.96 TeV 0-31, in fixed-target p- 
nucleus collisions with beam energies from 800 to 920 GeV [5[, and recently in pp collisions 
at 7 TeV at the Large Hadron Collider (LHC) ||6tiT3l. As the expected cross sections can be 
calculated in perturbative quantum chromodynamics (QCD), the comparison between data 
and predictions provides a critical test of next- to-leading-order (NLO) calculations |fl4llT5| . 

Considerable progress has been achieved in understanding heavy-quark production at Teva- 
tron energies, largely resolving earlier discrepancies in which theoretical predictions were sig- 
nificantly below observed production rates [TBI . However, substantial theoretical uncertain- 
ties on production cross sections remain due to the dependence on the renormalization and 
factorization scales. Measurements of b-hadron production at 7 TeV represent a test of the- 
oretical approaches that aim to describe heavy-flavor production at the new center-of-mass 
energy |[T6"IIT7|| . Furthermore, understanding the production rates for b hadrons represents an 
essential component in accurately estimating heavy-quark backgrounds for various searches, 
such as H° — > bb and supersymmetric or exotic new physics signatures with b quarks. 

This Letter presents the first measurement of the production cross section of a b baryon, A^, 
from fully reconstructed ]/tp A decays in pp collisions at a/s = 7 TeV and complements the mea- 
surements of B + [6J, B° [|7j, and Bg [9] production cross sections also performed by the Compact 
Muon Solenoid (CMS) experiment at the LHC [18]. The comparison of baryon production rela- 
tive to meson production resulting from the same initial b-quark momentum spectrum allows 
for tests of differences in the hadronization process. Such differences are particularly interest- 
ing in the context of heavy-baryon production in relativistic heavy-ion collisions, where the 
medium could significantly enhance the production of heavy baryons relative to mesons I1T9T - 
|2T| . Furthermore, the pp initial state at the LHC allows tests of baryon transport models, which 
predict rapidity-dependent antibaryon/baryon asymmetries, in contrast to baryon-antibaryon 
pair production, which typically results in equal yields |l22ll23| . Measurements of the A^ to A^ 
cross-section ratio, c(Ab)/0"(Ab), as functions of p Ab and \y Ab \ allow for the first test of such 
models with heavy-quark baryons at \/s = 7 TeV. 

Events with A^ baryons reconstructed from their decays to the final state ]/tp A, with J/tp — > 
and A — > pzr, are used to measure the differential cross sections da/dp Ah x B(A^ — > 
]/ipA),dtr/dy Ah x £>(Ab —> J/^A), and c"(Ab)/c"(Ab) with respect to the transverse momentum 
p Ab and the rapidity |y Ab |, as well as the integrated cross section times branching fraction for 
p Ab > 10 GeV and \y Ab \ < 2.0. The cross section times branching fraction is reported instead of 
the cross section itself because of the 54% uncertainty on £>(Ab — > J/?/ 7 A) |24"|. The cross section 
times branching fraction measurements are averaged over particle and antiparticle states, while 
the ratio is computed by distinguishing the two states via decays to p or p, respectively. 

2 Detector 

The data sample used in this analysis was collected by the CMS experiment in 2011 and cor- 
responds to an integrated luminosity of 1.86 ± 0.04 fb _1 |25| . A detailed description of the 
detector may be found elsewhere [18]. The main detector components used in this analysis are 
the silicon tracker and the muon detection systems. 

The silicon tracker measures charged particles within the pseudorapidity range < 2.5, 
where f] = — In [tan (9/2)] and 9 is the polar angle of the track relative to the counterclock- 
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wise beam direction. It consists of 1440 silicon pixel and 15 148 silicon strip detector modules 
and is located in the 3.8 T field of the superconducting solenoid. It provides an impact param- 
eter resolution of about 15 /mi and a pj resolution of about 1.5% for particles with transverse 
momenta up to 100 GeV. Muons are measured in the pseudorapidity range \rj\ < 2.4, with de- 
tection planes made using three technologies: drift tubes, cathode strip chambers, and resistive 
plate chambers. Events are recorded with a two-level trigger system. The first level is com- 
posed of custom hardware processors and uses information from the calorimeters and muon 
systems to select the most interesting events. The high-level trigger processor farm further 
decreases the event rate from about 100 kHz to around 350 Hz before data storage. 

3 Event selection 

Early data taking conditions in 2011 utilized a loose dimuon trigger with the following require- 
ments. Events are selected requiring two oppositely charged muons with dimuon transverse 
momentum greater than 6.9 GeV. Displaced muon pairs from long-lived b-hadron decays are 
preferentially selected by further requiring a transverse separation from the mean pp collision 
position ("beamspot") greater than three times its uncertainty, where the uncertainty incorpo- 
rates the vertex and beamspot measurements. Also required at the trigger level are a dimuon 
vertex fit confidence level larger than 0.5% and cos a > 0.9, where a is defined as the angle 
in the plane transverse to the beams between the dimuon momentum and the vector from 
the beamspot to the dimuon vertex. The dimuon invariant mass m fi + }i - is required to satisfy 
2.9 < m^+p- < 3.3 GeV. For the later 46% of the dataset, the trigger was tightened by increasing 
the dimuon vertex fit confidence level threshold to 10% and imposing kinematic requirements 
of pj > 3.5 GeV and \rf \ < 2.2 for each of the muons. The remaining 2011 data were recorded 
with even tighter triggers and are not used in the analysis. 

Muon candidates are fully reconstructed by combining information from the silicon tracker |26"| 
and muon detectors, and are required to be within the kinematic acceptance region of p\ > 
3.5 GeV and \rjl*\ < 2.2. Muon candidates are further required to have a track x 2 per degree 
of freedom <1.8, at least 11 silicon tracker hits, at least two hits in the pixel system, and to be 
matched to at least one track segment in the muon system. Multiple muon candidates are not 
allowed to share the same muon track segments [27]. 

Opposite-sign muon pairs are fit to a common vertex to form J/tp candidates, which are required 
to be within 150 MeV of the world-average ]/ip mass [24]. The ]/rp candidates are also required 
to have pj greater than 7 GeV, a dimuon vertex fit confidence level larger than 0.5%, cos a > 
0.95, and a transverse separation of the vertex from the beamspot greater than three times its 
uncertainty. 

The A candidates are formed by fitting oppositely charged tracks to a common vertex. Each 
track is required to have at least 6 hits in the silicon tracker, a x 2 P er degree of freedom <5, 
and a transverse impact parameter with respect to the beamspot greater than 0.5 times its un- 
certainty. The proton candidate, identified as the higher-momentum track, is required to have 
pj > 1.0 GeV. Misassignment of the correct proton track is found to be negligible from sim- 
ulation. The reconstructed A decay vertex must have a x 2 P er degree of freedom <7 and a 
transverse separation from the beamspot at least five times larger than its uncertainty. The in- 
variant mass m p7r is required to be within 8 MeV of the world-average A mass [24]. Candidates 
are rejected if m n + n - is within 20 MeV of the world-average Kg mass [24]. 

The At, candidates are formed by combining a ]/ip candidate with a A candidate. A vertex- 
constrained fit is performed with the two muons and the A candidate, with the invariant 
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masses of the J/ip and A candidates constrained to their world-average values [24 J. The Ab 
vertex fit confidence level is required to be greater than 1% and the reconstructed Ab mass 
must satisfy 5.2 < mj/^A < 6.0 GeV. Multiple Ab candidates are found in less than 1% of the 
events with at least one candidate passing all selection criteria. In those cases, only the can- 
didate with the highest Ab vertex fit confidence level is retained. The /Wj/^a distributions for 
selected Ab and Ab candidates are shown in Fig.jl] 
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Figure 1: Fit results for the niy^ A distribution for Ab (left) and Ab (right) for pj b > 10 GeV and 
\y Ab | < 2.0, where the dashed line shows the background fit function, the solid line shows the 
sum of signal and background, and the points indicate the data. 



4 Efficiency determination 

The efficiency for triggering on and reconstructing Ab baryons is computed with a combina- 
tion of techniques using the data and large samples of fully simulated Monte Carlo (MC) sig- 
nal events generated with pythia 6.422 |28], decayed by evtgen J29l, and simulated using 
Geant4 [30 1 . The efficiency is factorized according to 

fc trig fc trig fcreco treco t t rig fc sel' 

where each term is described below. The trigger (ejjL) an< ^ muon-reconstruction efficiencies 
(efeco) are obtained from a large sample of inclusive ]/tp — > decays in data using a "tag- 

and-probe" technique similar to that described in Ref. |3lT , where one muon is identified with 
stringent quality requirements and the second muon is identified using information either ex- 
clusively from the tracker (to measure the trigger and offline muon-identification efficiencies) 
or from the muon system (to measure the trigger and offline tracking efficiencies). While, in 
principle, the inclusive ]/ip — > sample can include signal events, which could bias the 

measurement, in practice the fraction is negligibly small and provides an unbiased measure- 
ment of the muon efficiencies. 

For the portion of the trigger efficiency that depends on single-muon requirements (e|r ig ), the 
efficiency for a given Ab event is computed as the product of the two single-muon efficien- 
cies. However, the trigger efficiencies for dimuon events where the muons bend toward each 
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other are up to 30% lower than for events where the muons bend away from each other for 
certain portions of the detector. This inefficiency arises when the muon trajectories cross in the 
muon system, and one of the candidates is rejected because of shared hits. To account for this 
effect, the trigger efficiencies for muons that bend toward and away from each other are com- 
puted separately in data and the appropriate efficiency is applied to each class of signal events. 
This procedure naturally accounts for the correlations between the two single-muon efficien- 
cies, as confirmed in simulation. The portions of the trigger efficiency that depend on dimuon 
quantities (e^L) are measured from an inclusive ]/ip sample collected with triggers where only 
single-muon requirements are applied. 

The probabilities for the muons to lie within the dimuon kinematic acceptance region (A) and 
for the At and Ab candidates to pass the selection requirements (e^) are determined from 

the simulated events. To minimize the effect of the pythia modeling of the pj h and |j/ Ab | 
distributions on the acceptance and efficiency calculations, the simulated events are reweighted 
to match the kinematic distributions observed in the data. The simulated events used for the 
efficiency calculations have also been reweighted to match the measured distribution of the 
number of pp interactions per event (pileup). On average, there are six pileup interactions in 
the data sample used in this analysis. The efficiencies for hadron track reconstruction 11321 , At, 
reconstruction 1133 1 , and fulfilling the vertex quality requirements are found to be consistent 
between data and simulation. 

The total efficiency of this selection, defined as the fraction of Ab — > ]/ipA with ]/xp — » p + p~ 
and A — > p7i decays produced with p Ab > 10 GeV and |y Ab | < 2.0 that pass all criteria, is 0.73%. 
The efficiency ranges from 0.3% for p Ab 10-13 GeV to 4.0% for pj h > 28 GeV, with the largest 
losses due to the A reconstruction (10-16% efficiency), the dimuon kinematic acceptance (12- 
63%), and the displaced dimuon trigger requirements (33-56%). The efficiencies in bins of pj h 
and |y Ab | are shown in Table [l] 

To measure the ratio of antiparticle to particle cross sections c"(Ab)/c(Ab), only the ratio of 
the Ab and Ab detection efficiencies is needed. Many of the efficiency contributions cancel in 
the ratio, including all the ]/ip and p. efficiencies since the particle and antiparticle states are 
indistinguishable. However, the A and A reconstruction efficiencies differ because of different 
interaction cross sections with the detector material; the p are more likely to suffer a nuclear 
interaction and be lost, resulting in an efficiency that is on average 13% lower for Ab than for 
Ab, as shown in Table [2] The ratio of the Ab and Ab selection efficiencies is calculated from 
simulation as described above for the combined sample, where the simulation modeling of the 
detector interactions is validated by comparing the number of hits reconstructed on tracks with 
that observed in data. The uncertainty on the amount of detector material and the appropri- 
ateness of simulated interaction cross sections are considered as systematic uncertainties, as 
described in Section [7] 

5 Fitting procedure 

The backgrounds are dominated by nonprompt ]/xp production from b hadrons. The dimuon 
invariant-mass distribution in data confirms that the contamination from events containing a 
misidentified ]/xp is negligible after all selection criteria have been applied. Background events 
are distinguished from signal by their reconstructed myy A distribution, which is found to be in 
good agreement between data away from the signal peak and simulated b — > ]/ipX events. The 
Ab proper decay length distribution in data confirms that the background events arise from 
long-lived b hadrons, and therefore offers no additional discriminating power between signal 
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and background. The measured m^ n distribution shows a purity of 77% genuine A events after 
applying the full selection criteria, while the m n + n - distribution confirms that more than 99.9% 
of the Kg background is rejected by the kaon mass-window veto. 

The Ab yields are extracted from unbinned extended maximum-likelihood fits to the ot^a 

distribution in bins of p Ah and |y Ab | defined in Table [lj In each bin, the signal is described 
by a double-Gaussian function with resolution parameters fixed to values found when fitting 
simulated signal events and means set to a common value left free in the fit. The background 
shape is modeled with a third-order polynomial, whose parameters are left free to float inde- 
pendently in each bin. The ratio of antiparticle to particle yields is obtained by simultaneously 
fitting the Ab and Ab mass distributions, with resolution parameters fixed from the fit to the 
combined Ab and Ab simulated sample and common mean allowed to float. The background 
shapes are fit with separate third-order polynomials, whose parameters are left free in the fit. 
The signal mass resolution varies as a function of \y Ab \, ranging from a mean of 11 MeV for 
central Ab to 27 MeV for forward Ab events. 



6 Results 

The fitted signal yields in each bin of p Ab and \y Ah | are summarized in Table [l] Figure [l] shows 
the fits to the niyy A distributions for Ab and Ab candidates in the inclusive sample with p Ab > 
lOGeV and \y Ab \ < 2.0. The total number of signal events extracted from an inclusive fit is 
1252 ± 42, where the uncertainty is statistical only. 



The Ab differential cross section times branching fraction is calculated in bins of p Ab as 

and similarly for \y Ab |, where n S i g is the fitted number of signal events in the given bin, e 
is the average efficiency for signal Ab and Ab baryons to pass all the selection criteria, £ is 
the integrated luminosity, Ap Ab is the bin size, and B is the product of branching fractions 
BQ/rp -4 = (5.93 ± 0.06) x HT 2 and B(A -» pn) = 0.639 ± 0.005 BSJ. The additional 

factor of two in the denominator accounts for our choice of quoting the cross section for Ab 
production only, while n S i g includes both Ab and Ab- The efficiencies are calculated separately 
for each bin, always considering only baryons produced with \y Ab | < 2.0 for p Ab bins and 
p Ab > lOGeV for \y Ab | bins, and taking into account bin-to-bin migrations (0-2%) because of 
the finite resolution on the measured p Ab and \y Ab |. Equal production of Ab and Ab is assumed 
for the efficiency, as predicted by pythia and as is consistent with our measurement. 

The measured differential cross sections times branching fraction versus p Ab and \y Ah \ are 
shown in Fig. [2] and Table [lj They are compared to predictions from the NLO MC generator 
POWHEG 1.0 with the hvq package 11341 1351 using a b-quark mass = 4.75 GeV, renormaliza- 

tion and factorization scales ]i = \Jnt^ + p\, CTEQ6M parton distribution functions [36], and 
PYTHIA 6.422 [28 1 for the parton hadronization. The uncertainty on the predicted cross sec- 
tion is calculated by varying the renormalization and factorization scales by factors of two and, 
independently, by ±0.25 GeV. The largest variation in each direction is taken as the uncer- 
tainty. The data are also compared to the PYTHIA 6.422 prediction, using a b-quark mass of 
4.80 GeV, CTEQ6L1 parton distribution functions, and the Z2 tune l|37| to simulate the under- 
lying event. No attempt has been made to quantify the uncertainty on the pythia predictions. 
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6 Results 



Table 1: A b + signal yield n S i g , efficiency e, and measured differential cross sections times 
branching fraction dc"/dpy b x £>(At, — > ]/ipA) and dcr/dy Ab x B(A b — > J/t/?A), compared to 
the POWHEG |34l [35] and PYTHIA |28l predictions. The uncertainties on the signal yields are 
statistical only while those on the efficiencies are systematic. The uncertainties in the measured 
cross sections are statistical and systematic, respectively excluding the common luminosity 
(2.2%) and branching fraction (1.3%) uncertainties. The POWHEG and PYTHIA predictions also 
have uncertainties of 54% due to B(A b — > ]/xpA), which are not shown. 
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The measured pj spectrum falls faster than predicted by POWHEG and PYTHIA, while the |y| 
spectrum shape is in agreement with the predictions within uncertainties, as illustrated in the 
data-to-POWHEG ratio plots shown in the lower panels of Fig. [2] The integrated cross section 
cr(pp -> A b X) x B(A b -4 ]/ipA) for pj b > 10 GeV and |y Ab | < 2.0, calculated as the sum 
over all pj bins, is 1.16 ± 0.06 ± 0.12 nb, where the first uncertainty is statistical, and the sec- 
ond is systematic. For the total cross section result, the highest pj b bin is fit without an upper 
bound and has a yield of 97.0 ± 13.2 events. The total cross section measurement is in good 
agreement with the prediction from PYTHIA of 1.19 ± 0.64nb and higher than the prediction 
from POWHEG of 0.63^37 nb, where the uncertainties are dominated by the 54% uncertainty 
on B{A h -> ]/tpA) [24]. ' 

This result can be compared to previous CMS measurements of B + 0, B° 0, and B° [Q pro- 
duction at = 7TeV. To facilitate the comparison, the B + and B° results are taken for the 
range p T > 10 GeV. Simulated events are generated with MC@NLO [38 1 with = 4.75 GeV 
and CTEQ6M parton distribution functions to determine the fraction of B + , B°, and B° events 
within the pj and \y B \ ranges used for their respective measurements with the pj > 10 GeV 
and |y| < 2.0 requirements used in this analysis. Scaling by the appropriate ratio and using 
the world-average values of B(A b -4 ]/tpA) = (5.7 ± 3.1) x 10~ 4 and B(B° S -4 J/iptp) = (1.4 ± 
0.5) x 10~ 3 ||24| , we determine the following cross sections for p T > 10 GeV and |y B | < 2.0: 
cr(pp -4 B+X) = 6.7 ± 1.0 ^b; a(pp -4 B°X) = 6.7±0.8pb; a{pp -4 B°X) = 2.5 ± 1.0 yb 
and cr(pp — > AbX) = 2.1 ± 1.1 pb, where the uncertainties are the quadrature sum of the 
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Figure 2: Upper: Measured differential cross sections times branching fraction da/ dp T h x 
B(A b — > ]/tpA) (left) and da/dy Ab x £>(A b —> J/t/?A) (right) compared to the theoretical predic- 
tions from PYTHIA and POWHEG. The inner error bars correspond to the statistical uncertainties 
and the outer ones represent the uncorrelated systematic uncertainties added in quadrature to 
the statistical uncertainties. The dashed lines show the uncertainties on the POWHEG predic- 
tions. Overall uncertainties of 2.2% for the luminosity and 1.3% for the ]/tp — > and 
A — > pre branching fractions for the data are not shown, nor is the 54% uncertainty due to 
£>(A b — > J/rpA) for the PYTHIA and POWHEG predictions. Lower: The ratio of the measured 
values to the POWHEG predictions. The error bars include the statistical and uncorrelated sys- 
tematic uncertainties on the data and the shape-only uncertainties on the POWHEG predictions. 



statistical and systematic components. No uncertainty has been included for the phase-space 
extrapolation based on MC@NLO [38J. The large systematic uncertainties for c(pp — > BgX) and 
C"(pp — > AbX) are dominated by the poorly known branching fractions B(A b — > ]/tpA) and 
£?(Bg — > ]/ip(p), respectively. The ratios among the four results are in good agreement with the 
world-average b-quark fragmentation results |24|. 

The world-average b-quark fragmentation results assume that the fractions are the same for b 
jets originating from Z decays at LEP and directly from pp collisions at the Tevatron. However, 
measurements of /A b performed at LEP ||59ll40l and at the Tevatron |j4H| show discrepancies. A 
recent result [42] from the LHCb Collaboration measures a strong pj dependence of the ratio of 
Ab production to B-meson production, /a,,/ {fu + fi), with f^ b = B(b — > A b ) and fq = Bib — > 
Bq ) . Larger f^ b values are observed at lower pj, which suggests that the discrepancy observed 
between the LEP and Tevatron data may be due to the lower pj of the A^, baryons produced at 
the Tevatron. 



A comparison of this and previous CMS results for b-hadron production versus pj is shown in 
the left plot of Fig. |3} where the data are fit to the Tsallis function [43J, 



1 ^ 
N dp T 



Cp T 



1 + 



Pj + m 2 



m 



nT 



(3) 



Here C is a normalization parameter, T and n are shape parameters, m is the mass of the b 
hadron and N is the b-hadron yield. The statistical and bin-to-bin systematic uncertainties are 
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used in the fits. The T parameter represents the inverse slope parameter of an exponential, 
which dominates at low pj. Since our data do not constrain that region well, T is fixed to the 
mean value found from fitting the B + and B° distributions, where the pj threshold is lowest. 
The result of T = 1.10 GeV is used to obtain the following values of the n parameter, which 
controls the power-law behavior at high p T : n(B+) = 5.5 ± 0.3, n(B°) = 5.8 ± 0.3, n(B°) = 
6.6 ± 0.4, and n(A\ } ) = 7.6 ± 0.4. The larger n value for Aj, indicates a more steeply falling pj 
distribution than observed for the mesons, also suggesting that the production of A^ baryons, 
relative to B mesons, varies as a function of pj, with a larger A^/B ratio at lower transverse 
momentum. The right plot of Fig. [3] shows the p A b spectrum shape compared to B + and B°, 
where the distributions are normalized to the common bin with pj = 10—13 GeV. 




b-hadron p t (GeV) b-hadron p t (GeV) 

Figure 3: Comparison of production rates for B + [6|, B° [7|, B° and A^ versus pj. The left 
plot shows the absolute comparison, where the inner error bars correspond to the total bin-to- 
bin uncertainties, while the outer error bars represent the total bin-to-bin and normalization 
uncertainties added in quadrature. Fits to the Tsallis function [43] for each distribution are 
also shown.The overall uncertainties for Bg and Ab are dominated by large uncertainties on 
£>(Bg — > ]/ip<p) and £>(Ab — > ]/tpA), respectively. The right plot shows a shape-only comparison 
where the data are normalized to the 10—13 GeV bin in pj and the error bars show the bin-to- 
bin uncertainties only. B^ is omitted because the 10—13 GeV bin is not available for the common 
normalization. 



The ratio c"(Ab)/V(Ab) is calculated in bins of p T b or \y Ab | as 

c(Ab)MAb) = ^| x ^1, (4) 

where n A ^ and n A ^ are the antiparticle and particle yields in a given bin, and e(Ab) and e(Ay ) ) 
are the particle and antiparticle efficiencies for a given bin, always considering only baryons 
produced with \y Ab \ < 2.0 for p Ab bins and p Ah > 10 GeV for |y Ab | bins. The results versus p Ab 
and \y Ah \ are shown in Fig.|i]and Table|2] The ratio <j(Ab)/c"(Ab) is found to be consistent with 
unity and constant as a function of both p Ab and \y Ab \, within the uncertainties, as predicted 
by POWHEG and pythia. Therefore, no evidence of increased baryon production at forward 
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pseudorapidities is observed within the available statistical precision for the kinematic regime 
investigated. The integrated (r(A b ) /cr(A b ) for pj h > lOGeV and \y Ab \ < 2.0 is 1.02 ± 0.07 ± 
0.09, where the first uncertainty is statistical and the second is systematic. 
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Figure 4: Measured c(Ab)/fj(Ab) (points) versus pj b (left) and |y Ab | (right), compared to the 
theoretical predictions from PYTHIA (red dashed line) and POWHEG (blue solid line). The inner 
error bars correspond to the statistical uncertainties, and the outer error bars represent the 
uncorrelated systematic uncertainties added in quadrature to the statistical uncertainties. The 
dashed blue lines show the uncertainties of the POWHEG predictions. 



7 Systematic uncertainties 

The cross section is affected by systematic uncertainties on the signal yields and efficiencies that 
are uncorrelated bin-to-bin and can affect the shapes of the distributions, and by the uncertain- 
ties on branching fractions and integrated luminosity, which are common to all bins and only 
affect the overall normalization. The uncertainties on the signal yields arise from the following 
sources: 

• Signal shape uncertainty (1-6%): evaluated from the variations when floating the 
means of the two Gaussians (set to a common value) in data or by using a single 
Gaussian shape. 

• Background shape uncertainty (1-2%): evaluated from the variation when using a 
second-order polynomial, exponential, or third-order polynomial fit in the restricted 
range 5.4-6.0 GeV. 

• Final-state radiation (0-1%): evaluated by removing it from the simulation and tak- 
ing half of the difference in the results. 

The uncertainties on the efficiencies arise from the following sources: 

• Pion/ proton/ A reconstruction efficiency uncertainty (8%): evaluated by varying the 
simulated detector material [44J, alignment, and beamspot position, and by varying 
the reconstruction cuts, by using different event simulations, and comparing the 
measured A lifetime (33|, which is sensitive to the efficiency correction, to the world- 
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7 Systematic uncertainties 



Table 2: Uncorrected signal yield ratio n^/n^, efficiency ratio e(A b )/e(A b ),, and efficiency- 
corrected ratio a(A h )/a(A b ), compared to the POWHEG EHISHl and PYTHIA |28| predictions 
for the corrected ratio. The uncertainties in the corrected ratio are statistical and systematic, 
respectively. The uncertainties on the uncorrected yield ratio are statistical only and on the 
efficiency ratio are systematic only. 
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Data 
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average value |24|. 

• Tag-and-probe statistical uncertainties (4-6%): evaluated by propagating statistical 
uncertainties from the data-driven determination of the single-muon efficiencies. 

• Tag-and-probe systematic uncertainties (1-7%): evaluated as the difference between 
the true efficiency in simulation and the efficiency calculated with the tag-and-probe 
procedure applied to simulated events. 

• Statistical precision of the simulated event samples (3-4%): calculated for the dimuon 
acceptance and reconstruction efficiencies. 

• Simulation modeling of the Ab kinematic distributions (0-5%): evaluated as half of 
the difference due to the kinematic re weighting. 

• Geant4 p cross section (1^%): evaluated by considering an alternative cross sec- 
tion model in GEANT4 (CHIPS) 11301 for p cross sections for interacting with material 
in the detector [45J and taking the difference in the efficiency as a systematic uncer- 
tainty. 

• Unknown Ab polarization (1-4%): evaluated by generating samples of events with 
the Ab spin fully aligned or anti-aligned with the normal to the plane defined by the 
Ab momentum and the pp beam direction in the laboratory frame and taking the 
average difference in the efficiency when compared to the nominal analysis, which 
is performed with unpolarized simulated events. 

• Pileup (0-4%): evaluated by varying the number of pileup interactions in simulated 
events by the uncertainty of the measured pileup interaction distribution. 

• Muon kinematics (0-2%): evaluated as the difference in the simulated efficiency 
when reweighting the muon pj to match the distribution measured with muons 
from the inclusive J/t/> sample used in the tag-and-probe measurements. 

• Effect of events migrating between pj and y bins due to resolution (0-1%): evaluated 
as half of the correction deduced from simulated events. 

The bin-to-bin systematic uncertainty is computed as the sum in quadrature of the individual 
uncertainties and is summarized in Table [I] In addition, there are normalization uncertainties 
of 2.2% from the luminosity measurement fl25| and of 1.3% from the J/i/> — > and A — > 

p7i branching fractions [24 j. For the total cross section result computed from the sum of pj 
bins, only the signal and background shapes, and the tag-and-probe and simulation statistical 
uncertainties are treated as uncorrelated. As bin-to-bin correlations cannot be ruled out for the 
remaining sources of systematic uncertainty, the contribution in each pj bin is added linearly 
to compute the sum to ensure that the uncertainty is not underestimated. 

Many of these systematic effects cancel in the c(Ab) / c(Ab) ratio measurement. The remaining 
uncertainties are from the signal shape (2-8%), background shape (1-3%), GEANT4 p cross 
section (1-7%), variation of detector material (5%), and statistical precision of the simulated 
samples (6-8%), which are evaluated as described above. The total systematic uncertainty is 
computed as the quadrature sum of the individual uncertainties and is summarized in Table [2] 

8 Conclusions 

In summary, the first measurements of the differential cross sections times branching fraction 
&a/dpj h x B(A b -> ]/ipA) and da/dy A * x B{A b -4 ]/tpA) for A b bary ons produced in pp col- 
lisions at y/s = 7TeV have been presented. The measurements are given for pj b > lOGeV and 
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\y Ab | < 2.0. The p T b distribution falls faster than both the measured pj spectra from b mesons 
and the predicted spectra from the NLO MC POWHEG and the leading-order MC pythia. The 
total cross section and rapidity distribution are consistent with both predictions within large 
uncertainties. The measured c(Ab) / c( Ab) ratio is consistent with unity and constant as a func- 
tion of both pj b and \y Ab \- 
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